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I. INTRODUCTION 



The quantum dissipation of a single mode harmonic oscillator coupled to a reservoir 
with an extremely large number of degrees of freedom has been intensely studied [l|-|5| 
in the past decades with various methods, such as master-equation approach, Langevin 
approach. By using the Markovian master-equation approach, the effect of dissipation on 
a macroscopic superposition of two coherent states of a harmonic oscillator was studied || . 
It shows that the superposition of the system will be transformed into a classical mixture 
state due to the effects of the environment 0, which is called decoherence process. The rate 
of the decoherence is proportional to the distance between the two states, which has been 
demonstrated in cavity QED experiment ||. The wave function structure of a single mode 
boson system plus a bath of many bosons is studied for both factorization case P,|TT)|] and 



partial factorization case pT| . It shows that when the Brownian motion effect is ignored 
in certain conditions the total wave function of the total system can be written in a form 
of a product of the bath and system components. Based on the result [O] , we have 



studied the decoherence process of excitons in an idealized quantum well placed in a lossy 



cavity [14]. The result shows that the coherence of excitonic superposition is reduced in an 



oscillating form due to the dissipation of the optical cavity. In references [|l"2"| , |i"3| , by solving 



the Heisenberg equations, the state vector of a driven-damped oscillator system evolved from 
an initial coherent state has also solved explicitly. 

As far as our knowledge is concerned, the dynamic evolution of a damped oscillator sys- 
tem is only investigated by implementing the Heisenberg-Schrodinger picture transformation 
(HSPT) [fnj-|I5||. Can we solve the Schrodinger equation directly and get the state vector of 
a damped harmonic oscillator system without HSPT? To answer this question, in this paper, 
we firstly apply normal ordering technique (NOT) to study the time evolution of a damped 
harmonic oscillator. The explicit form of the state vector of the total system evolved from an 
initially factorized coherent state is calculated. Besides, as one of the applications, we will 
devote ourselves to studying the decoherence behavior of the quantum dissipating system. 



For zero temperature case, the decoherence time of a superposition of two coherent states 
due to dissipation is given in a explicit form which agrees with the previous results P,|14 



The effect of finite temperature of the reservoir on the decoherence rate is also investigated. 

This paper is organized as following: in section II, we develop the normal ordering tech- 
nique to solve the dynamic evolution of a damped harmonic oscillator under certain initial 
condition. In section III, the decoherence of a superposition state of the single mode boson 
system is investigated. In subsection A, the decoherence process of the system in the case 
of zero temperature is studied for an factorized initial state. We find that, the decoherence 
occurs in a time scale which is sensitive to the distance of the initial superposed states. 
The temperature effect on the coherence evolution of a damped oscillator is investigated in 
subsection B by using NOT. 



II. NORMAL ORDERING TECHNIQUE 

The normal ordering technique is firstly introduced to study the dynamic evolution 
of a driven oscillator as well as that of two weakly coupled oscillators without dissipation. 
In this section we adopt the spirit of NOT to study the dissipation process of a single-mode 
boson system coupled with a reservoir. 

We consider a harmonic oscillator coupled to a reservoir which involves a large collection 
of systems with many degrees of freedom, that is a damped oscillator. The reservoir is 
modeled as a harmonic oscillator bath, which in practice can be the modes of the radiation 
or the quantized modes of elastic vibrations (phonons) in a solid. The Hamiltonian for the 
system plus bath can be described as 

H = frusta + h'^^Ujb'jbj + h s ^^{gjO i bj + H.c), (1) 

3 j 

where a(a') denotes the annihilation (creation) operator of the harmonic oscillator with 
frequency lu, and bj(bj) is the annihilation (creation) operator for the mode of the bath 
with frequency uij, which obey the boson commutation relations [6j , 6^] = 5 - gj denotes 
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coupling constant between the harmonic oscillator with frequency u and the oscillator mode 
of the reservoir with frequency Uj. The model described in eq.(l) can be exactly solved as 
long as the coupling constant gj and the spectrum density of the bath are specified explicitly 



IIJ . Here we adopt the spirit of the normal ordering method to study the dynamic evolution 
of a single mode oscillator with dissipation. 

The state vector of the whole system obeys the Schrodinger equation which has a solution 
of the form \if)(t)) = U(t) \if)(0)), where the time-evolution operator U(t) satisfies ihd t U{t) = 
HU(t), with the initial condition U(0) = 1. We assume that evolution operator has its 
normal order form U(t) = U^ n '{t). Since the normal form of any operator is unique, one can 
establish the one-to-one corresponding relationship between the normal ordered evolution 
operator U {n \t) and an ordinary function U (n \t), with U (n) (t) = (a, {(3j}\ U^(t) \{f3j} , a), 
where |{/3j}) = Ifo) ■ • ■ \Poo)- Such a corresponding relation defines a map , 

N- 1 : U«>(t) - U (n \t) = («, {&}| U^(t) , a) . (2) 

We may also define the inverse transformation K , 

K:F (n) (t) ^U {n \t) =U(t). (3) 

One can write the Schrodinger equation in the normal ordering form. With the operator 
one can further take diagonal coherent state matrix elements of both sides of the normal 

— (n) 

ordered Schrodinger equation and get a c-number equation of U , 



id t U {n) = [ua* + vfi U + JL) + ^ 9j 

3 i 3 

-*(a+4) + ? 9M (° + ^) ,I7( "' (4) 

Now we need to solve the equation for U (t) with the initial condition U (0) = 1. 
We let 

Z7 (n) = ex P [Aa*a + ]T Bfi]^ + ^ j ', V + E C ^ a + E D i a *Pi\> ( 5 ) 



Where the prime in denotes sum over index "j" and u j " with the condition j 7^ j'. Sub- 
stituting eq.(5) into eq.(4), we get the coupled equations for the time-dependent coefficients 

Bj = -iu^l + Bj), (6.a) 



Cj = -iujjCj - ig*{l + A), (6.b) 

i). = - iuDj _ ig .(i + _ ^ ' % , J B /ii , (6.c) 

3' 

A = -iu(l +A)—i 9i C 3i ( 6 - d ) 

3 

B^^-iujB^-igPj, (6.e) 

with the initial condition A(0) = £•■/(()) = Cj(0) = Dj(0) = 0. The solution of eqs.(6) is 
easily obtained by using Wigner-Weisskopf approximation UJ, 

A(t) = u(t) - 1, Bj = - 1, (7.a) 

C 3 {t)= Uj {t), D j (t)=v j (t), (7.b) 

where, for the bath with the frequency spectrum distribution p(uj) 

u (t) = e-^e^ 2 , (8.a) 

p i(u>j—ui)tp—jt/2 _ 1 

u (t) = -£* e -^* = — — , (8.b) 

i(u)j-!2)t p -yt/2 _ I 

„ ( t ) = -^--^ J , , ( 8 .c) 

cjj — uj + ry/2 



3,3 



here, 7 = 2np(u) \g(u>)\ 2 is the damping constant of the oscillator induced by the coupling 
to the environment, and uj = u + Au; is the renormalized physical frequency , with Au is 
the lamb frequency shift 



00 



PM \9(uj)\ 
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U) j — LU 



where "P" denotes the Cauthy principle part. 

If, as an example, the initial state of the total system is 

MO)) = |«> ® l{&}> , (10) 

where |{/%}) = Ylj denotes multimode coherent state of the bath. Then at any time t 
the total system will evolve into 



\m) 



3 



/ ^ fa-** + au 3 (t) + Y^v jtf (t)fy \), (11) 
where we have used the sum rule, 



which is derived from the normalized condition (ip(t)\ip(t)) = 1. We see from eq.(13), due 
to the bath fluctuation and the back-action of system on the bath, the state vector evolved 
from factorized initial state becomes fully entangled. If the Brownian effect caused by the 
terms Y2j fij v j{t) can b e ignored, the total state vector can be partially factorized fTI|| . We 



can further consider T = OK for the bath, that is, all the oscillator modes of the reservoir 
are in vacuum state initially. In cavity QED, this means that we neglect the background 
radiation of cavity. Then from eq.(13), the state vector can be simplified as 

\if)(t)) = \au(t))®\{au j {t)}), (13) 
and the corresponding sum rule is 

|«(t)| 2 + X>i(*)| 2 = l- (14) 
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III. DECOHERENCE OF SUPERPOSITION STATE OF THE HARMONIC 
OSCILLATOR DUE TO DISSIPATION 



Recently the mesoscopic superposition state of the cavity mode states M has been pre- 
pared experimentally. The phase shift between two components of the superposition can 
be well controlled by adjusting the interaction time between atoms and the cavity mode. 
Then the even and odd Schrodinger cats of cavity states were prepared in the experiment, 
with which a new scheme for logic qubit encoding are proposed to simplify the error 
correction circuits and improve the efficiency of the error correction in quantum compu- 
tation. Decoherence process of the superposed cavity state was also demonstrated in the 
experiment. Due to the presence of dissipation in the cavity, the coherent information of 
Schrodinger cat state will be lost within a time scale which is sensitive to the size of the 
meter states. 

As an application, in this section we will study the decoherence behavior of the damped 
oscillator described by eq.(l). Starting with eq.(ll), we will calculate the decoherence factor, 
which is defined as the coefficients of the off-diagonal element of the reduced density operator 
of system, for both T = OK and finite temperature of the heat bath. 

A. Zero-temperature case 

A single mode boson system coupled to a bath which is composed of a collection of 
many harmonic oscillators will result in the decohence process of the system, that is, a 
superposition state of the system will be transformed into a statistical mixture state due to 
the influence of the bath. If the initial state of the total system is 

|tf (0)> = (Cilai) + C 2 \a 2 }) ® |{0,}> , (15) 

that is, a superposition state for the system, and all the oscillator modes of the bath are 
in the vacuum states |{0j}) initially. Therefore from eq.(15), the state vector of the total 
system at time t can be written as 
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|*(t)) = dlamit)) I {a lUj (t)}) + C 2 |a 2 «(0) ® I {<W*)}>, (16) 

where the explicit form of w(t), are given in eq.(8), and eq.(ll.a), respectively. In 

addition, the sum rule 

$>,(t)| 2 = l-Kt)| 2 , (17) 

j 

has been used in eq.(16), which is derived from the condition, 

(*(t)|*(*)> = ICxI 2 + \C 2 \ 2 + dC^K) + C?C 2 < ai |a 2 ). (18) 

The dynamic decoherence process can be investigated quantitatively by calculating the 
reduced density matrix Tr jR (|\l/(t))(\l/(t)|) of the system at any time t. Then we will get the 
decoherence factor, which is defined as the coefficient of the non-diagonal elements of the 
reduced density matrix, such as 

F(t) = Y[(a 1 u j (t)\a 2 u j (t)) 

j 

= e (-5l°i| 2 -5l«2| 2 +< a 2)E J KW| 2 . ( 19 ) 

By using the sum rule of eq.(17), the decoherence factor becomes 

F(t) = e H |Ql|2 ^M 2 +^a2)(i-Ht)| 2 )_ ^20) 

The characteristic time of the decoherence of the superposition state is determined by 
the short time behavior of \F(t)\, that is, the norm of the decoherence factor. Within the 
time scale <C 1, the decoherence factor can be simplified as following, 

F(t) = e (-5l a il 2 -5l Q 2| 2 +^Q 2 )7^ (21) 

If we consider that ol\ = a, and a 2 = ae tA<p , where Aip is the phase shift of the initial 
superposed states. Then the characteristic time is determined as following, 

r" 1 = 2|a| 2 7 sin 2 (Aip/2) , (22) 



where \a\ 2 is the mean number of the oscillator. We can define the "distance" D = \a±— a 2 \ = 
2\a\ sin (Aip/2) between the two superposed states. Substituting D into eq.(22), we obtain 
that the characteristic time = jfi, where r p = I/7 is the life time of the oscillator due to 
its energy dissipation. Our result shows that the decoherence time is determined by both 
the phase difference of the meter states of initial superposition state and the mean number 
of the quanta of the single mode boson field for a fixed damping rate. For a special case 
Aip = n, which means that the system is prepared initially in an odd and even coherent 
states, the norm of decoherence factor of eq.(20) becomes 

\F(t)\ = e - 2 ^ 1 - e "' t \ (23) 

as a function of time. It shows that the coherence of the oscillator will decrease in the 
exponential decay rule. 



B. finite temperature case 

We assume that every oscillator mode of the bath is initially in thermal equilibrium state, 
and the single-mode boson is in a superposition of two coherent states. Then the density 
operator of total system at t — reads, 

p(0) = Ps <g> p b , (24) 

with 

Ps = (Ci|c*i) + C 2 |a 2 ))«a 2 | C* 2 + C?), (25a) 

and 

Pb = [ Y[d^,^—e-^ 2 /^ \(3 3 ) ((3,1 , (25.b) 
J j n \ n j) 

where (rij) denotes the mean occupation of the j-th oscillator mode with frequency Uj of 
the heat bath. From eq.(ll), the decoherence factor becomes 
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F(t) = e Hl«il a -sl<wl a +«I a a)(i-Kt)l a ) Y[ f jt (26) 

3 

with 

3 J n(n s ) 

By using the identity 

i y rf 2 /5exp(-A|/?| 2 + /i/5 + z/^) = ^exp(^), (28) 
with condition Re {A} > and arbitrary // and u, we obtain 

J = e -||aa-ai| 2 |M(*)| a |^(t)| a (nj>_ ( 2 g) 

Substituting of eq.(29) into eq.(26), we get 

,(-i|a 1 | 2 ~i| a2 | 2 +^a 2 )(l-|«(t)| 2 ) p --i| a2 -a 1 | 2 Kt)| 2 E,l^(t)| 2 (n J > 



p (-i| Q i| 2 -i|a2| 2 +^a 2 )(l-|«(t)| 2 ) p -i| a2 -«i| 2 Kt)|%(l- e -^) 



(30) 



where we have used the relation [12], 



^K-(t)| 2 (n,)=n(l-e-^), (31) 

here n = ( e fc s T — 1 ] , with is the Boltzmann constant and T is the absolute tempera- 
ture of heat bath. For low temperature, n ~ 0, the decoherence factor of eq.(30) will go back 
eq.(20). While at high temperature, n ~ as the previous subsection, the characteristic 
time Td of the decoherence of the superposition state is determined by calculating the norm 
of the decoherence factor within the time scale ft <C 1, that is, 

_P(t) = e (-il"i| 2 -f l a 2| 2 +^a 2 )7i e -i^l«2-ai|V_ (32) 

For oil = a, and 0:2 = ae tAlp , with A<£> is the phase shift between the "meter" states of the 
initial superposition state. Then the decoherence time is determined as 



1 = 2|«| 2 7 (l + ^)sin 2 (A^/2). (33) 
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As the previous subsection we substitute the "distance" D = \ai — a 2 \ = 2\a\ sin(A</?/2) be- 
tween the two superposed states into eq.(33) and obtain the decoherence time r d = — \^r- 
Comparing with the zero temperature case, we find that due to the finite temperature effect 
of bath the decoherence rate of the system becomes faster. 



IV. CONCLUSION 

The quantum dissipation process of a single-mode boson immersed in a bath of bosons is 
studied by using the normal ordering technique. The dynamic evolution of the total system 
is obtained and is used to study the decoherence behavior of the system both for the zero 
temperature case and for the finite temperature case. Due to the influence of the dissipation 
the coherence information of the system will be lost in a time scale which is dependent on 
the distance of the initial superposition state. For T ^ OK, the finite temperature effect of 
the bath will decrease the decoherence time. One can note that the normal ordering method 
can be also used to study the dynamic evolution of a driven-damped oscillator system. 
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